A direct experimental limit on neutron — mirror neutron oscillations 
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In case a mirror world with a copy of our ordinary particle spectrum would exist, the neutron n 
and its degenerate partner, the mirror neutron n', could potentially mix and undergo nn' oscillations. 
The interaction of an ordinary magnetic field with the ordinary neutron would lift the degeneracy 
between the mirror partners, diminish the n'-amplitude in the n-wavefunction and, thus, suppress 
its observability. We report an experimental comparison of ultracold neutron storage in a trap with 
and without superimposed magnetic field. No influence of the magnetic field is found and, assuming 
negligible mirror magnetic fields, a limit on the oscillation time Tnn' > 103 s (95% C.L.) is derived. 

PACS numbers; ll.30.Er, ll.30.Fs, 14.20.Dh, 28.20.-v 



The concept of a mirror world, as an attempt to restore 
global parity symmetry, has attracted interest since the 
1950's, started by the famous paper of Lee and Yang [l| 
and significantly expanded in the work of Kobzarev, 
Okun, and Pomeranchuk 0] . The mirror matter idea was 
first appUed to the Standard Model of particle physics 
in [3] ■ More recent overviews can be found in 0, [a| • The 
mirror world could hold a copy of the particle spectrum 
of our ordinary world. Matter and mirror matter would 
interact via gravity and present a viable explanation to 
the dark matter problem j^, 0, [1, H, • Besides grav- 
ity, other (new) interactions could show up in minute 
mixings of neutral matter particles — such as neutrinos, 
pions, kaons, or positronium (see TT] for e^e~) — and 
their degenerate mirror partners making oscillations be- 
tween them possible. Recently it was pointed out [l^ 
that no direct experimental limits exist on the oscilla- 
tion time Tnn' between ordinary matter neutrons (n) 
and the speculative mirror neutrons (n'). An indirect 



limit of the order Tn,,' > Is has been derived 



based on the search for neutron - antineutron (nn) os- 
cillations Fast nn' oscillations with Tnn' ~ 1 s, or 



at least much shorter than the neutron /3-decay lifetime, 
could explain 13, 3] the origin of ultra- high energy cos- 
mic rays above the Greisen-Zatsepin-Kuzmin (GZK) cut- 
off 1^ 17 1 . The viability of models and implications have 
been further discussed in [l^. 

Possible approaches to nn' oscillation experiments with 
sensitivities of several hundred seconds have been dis- 
cussed in [l^. One approach is to search for nn' oscil- 
lations by comparing the storage of ultracold neutrons 



(UCN) in vacuum in a trap in the presence and the 
absence, respectively, of a magnetic field. The essen- 
tial idea is that the neutron and mirror neutron states 
would be degenerate in the absence of a magnetic field 
and nn' transitions could occur. (The absence of mirror 
magnetic fields at the location of the experiment is as- 
sumed throughout this paper [131 ■) The interaction of the 
neutron with a magnetic field would lift the degeneracy 
and suppress the transition into a mirror neutron which, 
of course, does not interact with the ordinary magnetic 
field, nor with the trap via the ordinary strong interac- 
tion. Thus, the oscillation into mirror neutrons adds a 
loss channel for ultracold neutron storage. If nn' tran- 
sitions occurred, the storage time constant for ultracold 
neutrons in a trap with magnetic field would be longer 
than without magnetic field. One should note that this 
disappearance method only measures neutron loss as a 
function of applied magnetic field. A signal will not prove 
the oscillation into mirror neutrons, only that some mag- 
netic field dependent loss channel exists. By assuming 
that the nn' oscillation is responsible for the loss, limits 
can be set on Tnn'- One can imagine other exotic disap- 
pearance channels for the neutron, among which only the 
antineutron channel is tightly constrained ^14;] . 



The formulation of the nn' oscillation is analogous to 
the evolution of other simple two state systems such as 
spin i, K°K0, or nh mixing (see, e.g., [2l|, [13). The 
energy difference between neutron and mirror neutron 
states due to magnetic field interaction with the neutron 
magnetic moment a is ^lB. For convenience we adopt the 
notation of [13, and define a characteristic frequency 
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Lo = ^ which corresponds to half the energy splitting. 
The probability p for an UCN to be found as a mirror 
neutron after a time t can then be written as 



P{t) = 



Sin^(^l + {ujT^n'Y X t/Tnn') 



(1) 



The time t is limited by the free flight time tf between 
two wall collisions. The wall collision frequency is deter- 
mined by jj. The effective transition rate of UCN into 
mirror neutrons is then given by 



(2) 



For a real system the factors on the right hand side of 
Eq. ^ must be properly averaged over the distribution 
of flight times between collisions during the storage time 



Rs 



{p{tf))t 



(3) 



In experiments, one searches for a weak coupling, thus 
long T,in' : so in practice ujT^n' ^ 1 in Eq. ([T]) . Two 
limits are considered for Eq. ([1]): In the first case ("ti"), 
oJlitf ^ 1 (large B- field), many oscillations take place 
and the sin^(..) term of Eq. ([T|) is averaged to i because 
tf varies along UCN trajectories: 



^s,u - 



1 



(4) 



In the second case ("0"), ujotf <C 1 (small B-field), the 
n' component grows quadratically in time during the free 
flight: 



{tf)t 



(5) 



Also regular losses of UCN must be considered, such 
as absorption and upscattering (during wall interactions 
or in collisions with rest gas), trap leakage and /3 decay. 
All these loss mechanisms contribute to the UCN storage 
time constant Tgtoro of the system; Agtorc — 1 /''"store is the 
corresponding loss rate. Generally, the loss rate depends 
on UCN energy and for a spectrum of stored UCN the 
decay curve is a sum of exponentials. The total effect 
can be modeled by the relative populations Ci of differ- 
ent velocity classes, each with its own storage loss rate 
A'-*'. After storing an initial number n{t=0) of UCN for 
some time in a given magnetic field one will detect the 
number of surviving UCN 

n{ts) ^n{t^0)xJ2 exp[-(AijL + Rs) ts] ■ (6) 

i 

with the simple normalization condition Ci = 1. For 
measurements in the limits ti and (only the magnetic 



field is changed), the ratio of detected UCN becomes in- 
dependent of all the regular UCN loss mechanisms 



o/n 



exp[(i?sjX - Rs,a) ts 



(7) 



In the absence of other effects, neutron - mirror neutron 



oscillations lead to Ni 



o/n 



< 1. 



We have performed UCN storage experiments at the 
Institut Laue-Langevin using the experimental setup of 



the neutron EDM experiment [2J, [25|. A typical mea- 
surement cycle consists of (i) a filling time of 40 s with 
the beam switch connecting the storage chamber to the 
ILL PF2 EDM beam line ^ allowing unpolarized UCN 
to enter the storage volume, (ii) different storage times 
ts when the UCN isolation shutter to the storage cham- 
ber was closed and (iii) 40 s counting time with the UCN 
shutter open and the beam switch connecting the storage 
chamber to the ^He filled UCN detector [23|. The pres- 
sure inside the storage chamber was always better than 
10"'^ mbar in order to make sure that the nn' degeneracy 
is not lifted by the interaction of UCN with the rest gas. 

The UCN storage chamber has a volume V ^ 21 I and 
a surface area of A ~ 5400 cm^. The limit for stored 
UCN velocity is 4.1 m/s. From kinetic gas theory, the 
mean free path of UCN between wall collisions is ^ « 
0.16 m, the mean velocity is about 3 m/s [2^ and, thus, 
(tf) fa 0.053 s. One obtains lo ■ {tf) ^ 1 at a magnetic 
field of i? ~ 0.42 jiY; the limiting cases are obtained for 
magnetic fields of more than a few /xT (tJ,) and of less 
than a hundred nT (0), respectively. 

Different magnetic field configurations were used: up 
(B|), off (So = 0) and down {Bi). The strength of 
the magnetic field was adjusted by the current through 
the main magnetic field coil. The relevant measurements 
were taken at l-B^xl ~ 6/zT (100mA). The magnetic field 
{Bz^ along the main magnetic field direction) as a func- 
tion of the applied current was measured using the Hg 
cohabiting magnetometer [2^ for fields between 0.3 /j,T 
and 7 //T. For lower fields the Hg magnetometer could 
not be used. The Hg data set shows a perfectly linear 
dependence of the field on applied current and results in 
\Bzfi\ — 2.3 ± 2.6 nT when extrapolated to zero current. 
This value for the magnetic field along the main field axis 
indicates a residual absolute B-field below 13 nT because 
of the absence of a preferred spatial direction. The zero 
field Eq for the actual measurements was set by switching 
off the coil current and demagnetizing the four-layer mu- 
metal shield surrounding the storage chamber. We also 
used 3-axis fluxgate sensors directly above the storage 
chamber in order to verify that the residual B-field was 
sufficiently small for the purpose of our experiment. The 
B-field configurations of the experiment were very well 
reproducible, and in particular Bq within less than 1 nT. 
The direct limit on |i3o| obtained from the fluxgates is, 
however, somewhat weaker: it was found that the con- 
nectors of the commercial devices are slightly magnetic. 
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ts [s] 

[s] 


50 (a) 
73 ± 3 (a) 


50 (b) 
73 ± 3 (b) 


100 (a) 
123 ± 3 (a) 


175 (a) 
198 ± 3 (a) 


n{Bo) 


44317 ± 40 
44197 ± 53 
44128 ± 53 
44163 ± 38 


44363 ± 53 
44443 ± 53 
44316 ± 46 
44371 ± 35 


28635 ± 21 
28671 ± 30 
28596 ± 30 
28633 ± 22 


17015 ± 22 
17047 ± 31 
16974 ± 31 
17011 ± 22 




1.0035(13) 
1.001 


0.9998(15) 
9(10) 


1.0001(11) 


1.0002(18) 



TABLE I: Measured total UCN counts n normalized per cycle 
for the day sequences ((a), (b), see text) at different storage 
times ts (with systematic error) and magnetic field configu- 
rations. n{Bii) is the weighted average of n{B^) and n{Bi), 
and iVo/ti = n{Bo)/n{Bfi). 



leading to offset fields at the location of the sensor on 
the order of 25 nT. Although the residual field is proba- 
bly on the level of a few nT, we give a conservative limit 
of l-Bol < 50 nT, which is sufficient for our purpose here, 
i.e., for the limiting case "0". 

Most of the measurements were performed repeat- 
ing the sequence (a) {Bq, B-^, Bi, Bq, Bq, Bi, Bq, B-^) with 
field changes typically every 1.5 h during day time. The 
demagnetization procedure before Bq measurements took 
about half an hour. For a given B-field configuration 16 
UCN cycles were measured: 4 for each storage time of 
tg = 100 s, 50 s, 175 s, and again 100 s. Night runs were 
taken for longer periods at one B-field configuration with 
ts = 100 s. They were used to check on the long term 
stability of the system. It was found that drifts of the 
count rates were slow and on a level below 0.3% over 
several hours. This agrees in magnitude with changes in 
the reactor power, but a direct correlation could not be 
established. The count rate drifts were sufficiently slow 
to be averaged out in the day runs with frequent changes 
of the B-field configuration. For ts — 50 s some data 
was taken using another sequence (b), {B^ , Bi, Bq) while 
checking on an unexpected count ratio A^o/Ti (^^^ below). 

The time constant for UCN to leave the storage cham- 
ber with the shutter open was measured to be Top = 
11.4 ± 0.6 s. Mirror neutrons can leave the system also 
during filling and counting. The relevant average times 

in our storage system are thus longer than the times 
between closing and opening the UCN shutter. One can 
replace ts in Eq. ([7]) by and, because Rs changes only 
rather weakly with ts and Top is much smaller than ts, 
one finds to very good approximation t* — ts + 2 Top. We 
assign a conservative systematic error of ±3s to i*. 

For data analysis, each B-field configuration was first 
treated separately. The distributions of single-cycle 
counts n were found to be consistent with Gaussians with 
standard deviations a ~ with no additional systemat- 
ics. The counts per cycle for each configuration were thus 
statistically averaged, see Table HI The averaged n data 



show a presently unexplained tendency to a linear depen- 
dence on the magnetic field. The effect we wish to limit 
depends on \B\'^ (via w^), so the direct average of the -f 6 
and —6 /iT measurement values cancels the linear sys- 
tematic effect, leaving only the possible oscillation effect 
and any remaining quadratic systematic contributions. 
The averaged n data is then used to calculate the count 
ratios Nq/^^ according to Eq. fT]). For 50 s storage time 
and sequence (a), the count ratio Aq/i^ is larger than 1 by 
2.7 standard deviations, which led us to remeasure at this 
storage time (using sequence (b)) and corroborate that 
this deviation was a statistical fluctuation. Both ratios 
A'o/||(50 s) are given in the table along with the obtained 
average. The individual results for Nq^^^ show no signal 
within their respective sensitivities and, as they are inde- 
pendent, can be used in a combined analysis. Following 
Eq. ([7]), we write Aq/i^ — cxp[at<;] with a fit parameter 
a. The fit gives a = (5.38 ±5.78) x 10"^ s^^. We use this 
value at the limit of the experimentally measured range 
(t* = 198 s) to set the constraint on the neutrons which 
may have been lost, yielding 

7Vo/|x(t* = 198s) = 1.00106 ±0.00114. (8) 

Results with Nq/^^ > 1 are unphysical for nn' oscilla- 
tions. In order to derive a limit on Nq/-^^ we adopt the 
Bayesian approach described by the Particle Data Group 
(page 305 of Ref. [131 ; probability distribution set to zero 
for A"q/|| > 1) and obtain 

No/^i it*s = 198 s) > 0.99835 s (95%C.L.). (9) 

In order to derive the limit on Tnn/, the flight time 
distribution averages (t/)t, and (i^)*, are needed as ad- 
ditional input. A better determination than the one from 
the kinetic gas theory argument given above was obtained 
by Monte Carlo calculations using Geant4UCN [3l[ . The 
parameters of the simulation (mainly material properties, 
such as Fermi pseudo-potential, loss probability per wall 
reflection, and fraction of diffuse to specular reflection) 
have been tuned to reproduce measurements of the UCN 
beam energy spectrum and fiUing, storage, and empty- 
ing time curves of the apparatus. Excellent agreement 
with the observables is obtained, which justifies the ex- 
traction of the required flight time distributions. The 
averages are given in Table [ill The assigned systematic 
uncertainties were derived by varying the material pa- 
rameters of the simulation, the largest influence coming 
from the loss probabilities per wall collision. 

The limit on the oscillation time is obtained using the 
limit on lost neutrons, Eq. and the average free flight 
time values (at ts = 175 s, from Table HI)) in Eq. ([7]) and 
solving for r„n'. The systematic uncertainties are taken 
into account for the limit on T^n' by adding (or subtract- 
ing) them simultaneously in order for them to weaken 
the limit, i.e., {tf)t, = 0.0548s, {tj)t, = 0.00515s2 and 
t* — 195 s. With a reminder of the assumptions (negligi- 
ble mirror magnetic field (20| , no conventional strong or 
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ts [s] 


50 


100 


175 


{tf)ts [S] 


0.0498(5) 
0.00420(8) 


0.0515(5) 
0.00450(9) 


0.0543(5) 
0.00505(10) 



Polish Ministry of Science and Higher Education, con- 
tract No. 336/P03/2005/28, and the Swiss National Sci- 
ence Foundation #200020-111958. 



TABLE II: Results for {tf)t, and {t'j)t, using Monte Carlo 
distributions of flight times between wall collisions. Variation 
of parameters in the simulation is used to assign systematic 
uncertainties. 




90 100 110 120 130 

Oscillation Time t , [s] 



FIG. 1: (color online) The count ratio A^o/u ^ ^ function of 
the oscillation time (see Eq. (17|). The dashed region indicates 
the allowed region of Eq. (|9]). The error band of the curve dis- 
plays the systematic uncertainties, see text. The X indicates 
the point at which the limit on Tj^^i is evaluated. 



electro-magnetic interactions of the mirror neutron n', 
and degeneracy of n and n' in the gravitational field) we 
obtain the final result 



> 103s(95%C.L.). 



(10) 



Figure [T] displays the dependence of Nq/^^ on Tnn' (see 
Eq. ([7])) with the band indicating the influence of the 
systematic uncertainties. The 95% confidence limit of 
Eq. ^ is shown and the cross marks the point which 
determines the limit on the oscillation time. 

The result impacts the role nn' oscillations can play in 
the transport of ultra-high energy cosmic rays over large 



distances [1^ |15|, although, it may not completely rule 
out the nn'-explanation for events above the GZK cutoff. 
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